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Summary 


We  model  45  P  and  PP  waveforms  at  epicentral  distances  of  14°  to  40°,  primarily  sampling 
western  and  central  China,  to  determine  regionalized  P  wave  velocity  variations  for  the  crust  and 
upper  mantle  in  this  region.  The  association  between  these  variations  and  subsurface  tectonic 
features  is  established.  The  regionalized  velocity  models  are  used  to  generate  synthetic 
seismograms  for  explosion  and  earthquake  sources  at  regional  distances  to  examine  the  effect  of 
crust  and  mantle  lid  variations  on  the  propagation  of  the  Pn  and  Lg  phases. 

Observed  P  and  PP  waveforms  are  matched  by  one  dimensional  forward  modeling  using  the 
reflectivity  technique.  Our  approach  is  to  first  find  a  homogeneous  average  velocity  model  for 
China  which  can  match  the  observed  broadband  waveforms  ^tered  with  a  long-period  instrument 
response.  This  model  is  then  used  as  a  starting  model  for  the  matching  of  the  broadband 
waveforms  to  interrogate  the  lateral  velocity  variations.  Our  preferred  average  model  WCH  has  a 
50  km  thick  crust  with  a  velocity  of  6.4  km/s,  a  Pn  velocity  of  8. 15  km/s  and  a  low  velocity  zone 
between  100  km  and  175  km  depth.  It  has  intermediate  characteristics  between  an  active  tectonic 
region  and  a  stable  shield,  compatible  with  previous  S  wave  modeling  results.  Modeling  of  the 
broadband  waveforms  indicates  significant  lateral  deviations  from  model  WCH.  Paths  traversing 
eastern  China  require  a  constant  velocity  lid  of  about  100  km  thickness,  a  Pn  velocity  from  8.12- 
8.14  km/s  and  a  crustal  thickness  between  27-35  km.  For  northern  China  a  crustal  thickness 
between  43-50  km,  a  Pn  velocity  of  about  8.0  km/s  and  a  100  km  thick  lid  with  a  positive  velocity 
gradient  are  preferred.  For  the  Tibetan  plateau  we  observe  large  crustal  thickness  (50-60  km)  and 
Pn  velocity  (8.15-8.25  km/s)  in  combination  with  a  more  pronounced  low  velocity  zone  than  that 
of  the  average  model  WCH.  Lower  velocities  may  extend  below  200  km  depth.  Our  modeling 
indicates  a  ^in  lid  of  about  50  km  thickness  for  this  plateau.  We  infer  from  our  results  and 
previous  work  that  crustal  shortening  and  thickening  is  the  most  likely  process  responsible  for  the 
Tibetan  plateau’s  thick  crust  and  high  average  elevation.  We  suggest  that  the  Indian  plate  has 
probably  only  underplated  the  southernmost  part  of  the  plateau.  Mantle  convection  associated  with 
the  crustal  shortening  process  has  led  to  strong  lateral  heterogeneity  in  the  upper  mantle  under 
Tibet.  Higher  Pn  velocities  in  western  Tibet  might  indicate  that  this  region  is  presently  undergoing 
crustal  shortening. 

The  regiond  variations  in  lid  structure  predict  large  variability  in  regional  wave  amplitudes, 
particularly  for  the  Pn/Lg  discriminant.  This  ratio  may  vary  by  as  much  as  an  order  of  magnitude 
for  crustal  and  upper  mantle  velocity  variations  typical  for  China.  Short-period  Lg  is  found  to  have 
very  stable  propagation  characteristics;  high  frequency  Pn  rms  amplitudes  exhibit  a  more 
complicated  behavior.  For  any  of  the  regionalized  Chinese  velocity  models,  short-period  Pn/Lg 
ratios  for  the  explosive  source  are  consistently  higher  than  those  for  the  earthquake  source. 
However  differences  in  this  ratio  between  the  explosion  source  and  the  earthquake  source  are 
considerably  smaller  than  variations  in  Pn/Lg  ratio  resulting  from  propagation  differences. 
Discrimination  attempts  based  on  this  ratio  should  therefore  be  undertaken  with  care  and 
propagation  paths  for  both  sources  have  to  be  similar.  Broadband  Pn/Lg  ratios  for  the  explosion 
source  can  be  both  higher  or  lower  than  those  for  the  earthquake  source  making  this  passband 
poorly  suited  for  discrimination  purposes. 
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1.  Introduction 


Determination  of  upper  mantle  velocity  variations  in  Eurasia  has  long  been  problematic  due  to 
the  lack  of  seismic  stations  or  data  availability  in  China  and  the  former  Soviet  Union.  Early  P  wave 
velocity  studies  (e.g.  King  &  Calcagnile,  1976;  Given  &  Helmberger,  1980;  Vinnik  & 
Saipbekova,  1984)  were  limited  to  northwestern  Eurasia,  using  Russian  nuclear  explosion  data 
recorded  at  stations  in  western  Europe.  Recently,  preliminary  analyses  using  data  recorded  within 
the  C.I.S.  have  been  conducted  (Goldstein  et  al.,  1992;  Gamero  et  al.,  1992),  refining  earlier 
models  for  the  P  wave  velocity  structure  under  Russia.  Because  large  parts  of  Eurasia  are  devoid 
of  earthquake  sources  and  most  stations  have  been  confined  to  the  margins  of  this  continent, 
studies  of  central  and  east  Asia  have  mainly  used  multi-bounce  S  waves  or  surface  waves  which 
span  longer  distances  than  direct  P  and  S  waves.  Grand  &  Helmberger  (1985)  used  differential 
travel-times  and  general  characteristics  of  long-period  S,  SS  and  SSS  waveforms  from  southeast 
Asian  earthquakes  recorded  at  stations  in  western  Europe,  to  explore  radial  and  lateral  variations  in 
shear  wave  velocity  structure  in  Eurasia.  Lyon-Caen  (1986)  used  S  and  SS  waveforms  sampling 
the  Indian  Shield  and  the  Tibetan  plateau  to  determine  S  wave  velocity  variations  in  this  region. 
This  work  was  extended  by  Zhao  et  al.  (1991)  using  S,  SS,  Love  and  Pnl  waveforms  for  paths 
traversing  Tibet  and  eastern  China.  Many  surface  wave  studies  have  been  conducted  to  place 
constraints  on  the  shear  wave  velocity  structure  beneath  Eurasia,  particularly  for  the  Tibetan  plateau 
(e.g.  Chen  &  Molnar,  1975;  Chun  &  Yoshi,  1977;  Patton,  1980;  Romanowicz,  1982,1984;  Jobert 
et  al.,  1985;  Brandon  &  Romanowicz,  1986;  Bouijot  &  Romanowicz,  1992). 

Most  waveform  modeling  studies  assume  laterally  homogeneous  material  properties  to  justify 
a  one  dimensional  modeling  approach.  For  long-period  waveform  analysis  this  assumption  is  often 
met  and  one  dimensional  earth  models  have  given  satisfactory  fits  to  observed  waveforms  (e.g. 
Burdick,  1981;  Grand  &  Helmberger,  1984,  1985;  Lyon-Caen,  1986,  Lefevre  &  Helmberger, 
1989).  For  this  study  we  will  rely  on  the  same  assumption  on  a  path  by  path  basis;  however,  since 
broadband  data  are  much  more  sensitive  to  lateral  inhomogeneity  (Paulssen,  1987)  and  because  the 
tectonic  setting  in  southeastern  Eurasia  is  very  heterogeneous,  we  are  required  to  use  shorter  and 
more  geographically  restricted  paths  than  most  previous  studies.  This  is  only  viable  using  data 
collected  within  the  region  under  study. 

In  1986,  the  Chinese  Digital  Seismic  Network  (CDSN)  of  broadband  stations  was  deployed, 
which  together  with  the  relatively  high  seismicity  in  this  region  offers  for  the  first  time  an 
opportunity  to  study  broadband  P  and  PP  waveforms  at  close  enough  distances  to  justify  localized 
one  dimensional  modeling  efforts.  During  the  last  7  years  of  operation,  the  CDSN  stations  have 
accumulated  a  moderate  number  of  high  quality  P  waveforms  sampling  the  upper  mantle  under 
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China,  enabling  a  study  of  the  mantle  velocity  structure  with  relatively  high  density  path  coverage. 
In  late  1989  the  broadband  IRIS/IDA  network  began  deployment  in  the  C.I.S..  A  preliminary 
analysis  using  data  recorded  by  this  network  has  been  conducted  by  Goldstein  et  al.  (1992). 
However,  the  amount  of  high  quality  data  available  for  the  IRIS/IDA  stations  is  not  yet  large 
enough  for  a  thorough  study  of  the  upper  mantle  under  central  Asia.  We  will  therefore  concentrate 
on  resolving  lateral  variations  in  the  upper  mantle  P  wave  velocity  structure  under  China  using 
broadband  P  and  PP  waveform:  recorded  at  CDSN  stations  for  earthquake  sources  with  well 
constrained  moment  tensor  solutions. 

Analysis  of  the  velocity  structure  under  China  is  particularly  interesting  because  of  its  strong 
tectonic  heterogeneity  related  to  the  collision  of  the  Indian  and  Eurasian  plates  40-60  Ma  ago, 
which  initiated  the  Himalayan  orogeny.  The  Tibetan  plateau,  a  region  of  numerous  fold  belts  in 
southwest  China  (Figure  1)  with  an  average  elevation  of  5000  m  above  sea  level,  is  certainly  the 
most  impressive  feature  related  to  this  collision.  Several  deep  refraction  studies  (Yuan  et  al.,  1986; 
Mooney,  1992)  as  well  as  surface  wave  studies  have  indicated  an  average  crustal  thickness  of  70 
km  for  this  plateau.  The  fold  system  continues  into  Nepal  and  northern  India  and  is  bordered  to  the 
south  by  the  Indian  stable  shield.  The  northern  part  of  China  and  Mongolia  is  also  characterized  by 
several  fold  belts.  These  are  bordered  to  the  north  by  the  Siberian  platform  and  are  separated  from 
the  southwestern  fold  system  by  the  Tarim  basin  to  the  west  and  the  Sino-Korean  platform  to  the 
east  (Figure  1).  Crustal  thicknesses  range  from  50  km  in  northwest  China  to  35  km  for  the  eastern 
platforms  (e.g.  Mooney,  1992). 

Two  very  different  mechanisms  have  been  suggested  to  explain  the  evolution  of  the  Tibetan 
plateau.  Argand  (1924)  originally  proposed  that  the  whole  of  Tibet  is  underlain  by  the  northern  part 
of  the  Indian  plate.  Studies  supporting  this  theory  indicated  that  this  underthrusting  might  extend  as 
much  as  1000  km  into  Eurasia  (Barazangi  &  Ni,  1982;  Ni  &  Barazangi,  1983).  The  buoyancy  of 
the  thick  crust  thus  created  would  be  able  to  support  the  weight  of  the  overlying  mountain  range 
and  give  the  Tibetan  plateau  its  high  average  elevation.  In  contrast,  others  (e.g.  England  and 
Houseman,  1986)  suggested  that  Tibet's  crust  was  thickened  and  elevated  due  to  folding  and  thrust 
faulting  under  the  north-south  compressional  regime  resulting  from  the  collision.  Molnar  (1988) 
gave  an  extensive  review  of  geophysical  constraints  used  to  infer  the  evolution  of  the  Tibetan 
plateau  and  its  surrounding  regions.  Seismological  studies  of  the  lid  velocity  structure  under  Tibet 
play  an  important  role  in  resolving  this  evolution.  A  thick  high  velocity  lid,  associated  with  colder 
material,  is  viewed  as  evidence  for  underthrusting,  while  a  thin  or  non-existing  lid  points  to  the 
absence  of  an  underthrusted  lithospheric  plate  and  possibly  upwelling  warmer  mantle  material. 

Studies  of  P  wave  velocities  in  the  mantle  directly  beneath  the  Moho  have  found  Pn  values 
between  8.1  km/s  (Molnar  &  Chen,  1984)  and  8.43  km/s  (Barazangzi  &  Ni,  1982)  for  the  Tibetan 
plateau.  Holt  &  Wallace  (1990)  put  narrower  constraints  on  Pn  velocities  for  this  region;  8.1  -  8.25 
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km/s.  They  find  Pn  velocities  for  northern  India  and  southeastern  China  to  be  around  8.0  km/s. 
Holt  &  Wallace  (1990)  explain  the  sharp  velocity  jump  of  0.2  km/s  they  observe  for  the  transition 
zone  between  the  Indian  shield  and  the  southern  Tibetan  plateau  by  increased  pressures  in  the 
mantle  at  the  base  of  the  crust  due  to  the  doubling  of  the  crustal  thickness,  assuming  an 
unperturbed  temperature  regime.  For  southern  Tibet  they  find  that  a  100  km  thick  lid  with  a 
positive  lid  velocity  gradient  best  explains  their  data.  Holt  &  Wallace  (1990)  infer  that  the  southern 
half  of  Tibet  is  being  underthrust  by  the  Indian  lithosphere.  Lyon  Caen  (1986)  finds  a  thin  shear 
wave  velocity  lid  underlain  by  a  distinct  low  velocity  zone  for  the  central  and  northern  Tibetan 
plateau.  This  indicates  that  this  underthrusting  has  not  taken  place  under  the  whole  of  Tibet.  In 
contrast,  the  velocity  model  of  Zhao  et  al.  (1991)  for  similar  paths  through  Tibet  has  a  slower  lid 
than  Lyon-Caen's  (1986)  and  does  not  exhibit  a  pronounced  low  velocity  zone.  The  shear  wave 
velocity  model  of  Zhao  et  al.  (1991)  seems  to  be  at  least  partly  influenced  by  propagation  through 
the  Chang  Thang  province  in  north-central  Tibet  (Figure  1),  a  region  for  which  Bourjot  & 
Romanowicz  (1992)  indicated  a  thin  high  velocity  crust  together  with  a  low  Sn  velocity.  Ni  & 
Barazangi  (1983)  also  found  a  low  Q  value  for  this  region,  observations  that  can  be  correlated  with 
recent  volcanism  (Gansser,  1980). 

Less  is  known  about  upper  mantle  velocities  and  related  subsurface  tectonic  features  in  other 
parts  of  China.  Preliminary  S  wave  results  of  Zhao  et  al.  (1991)  for  eastern  China  indicate  a 
distinct  low  velocity  zone  in  the  lid  together  with  a  higher  Sn  velocity  than  for  their  Tibetan  model. 
The  focus  of  this  study  is  to  analyze  large  scale  lateral  variations  in  the  upper-mantle  P  wave 
velocity  under  China,  establishing  association  with  tectonic  features  related  to  the  Indian-Eurasian 
collision. 

The  1992  megaton  explosion  at  the  Lop  Nor  test  site  in  the  Tarim  basin  (Figure  1) 
emphasized  that  China  remains  an  important  region  for  nuclear  test  monitoring.  Pn  and  Lg  are 
often  used  as  underground  nuclear  explosion  discriminants  and  yield  estimators.  Variations  in 
crustal  and  upper  mantle  velocity  structure  have  pronounced  effects  on  the  propagation  of  these 
regional  phases.  For  this  purpose  it  is  important  to  understand  how  the  propagation  of  the  Pn  and 
Lg  phases  is  affected  by  variations  in  the  Chinese  crust  and  upper  mantle  velocity  structure.  This 
will  enable  discrimination  studies  to  correct  for  path  effects  and  thus  facilitate  discrimination  and 
improve  yield  estimates.  The  thickness  and  velocity  structure  of  the  crust  strongly  influence  the 
propagation  of  Lg.  For  explosion  sources,  the  excitation  of  Lg  is  governed  by  the  P  wave  velocity 
near  the  source  and  S  wave  velocities  in  the  uppermost  mantle  (Frankel,  1989).  The  ratio  of  these 
two  velocities  determines  how  much  energy  can  be  trapped  in  the  crust  (Xie  &.  Lay,  1993).  Along 
with  crustal  structure,  the  behavior  of  the  Pn  phase,  which  travels  just  below  the  Moho  also 
depends  on  lid  velocity  gradients,  with  less  Pn  energy  being  lost  by  radiation  into  the  mantle  for 
positive  velocity  gradients.  We  will  generate  synthetic  seismograms  at  regional  distances  to  assess 
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how  the  Pn  and  Lg  phases  are  influenced  by  our  models  of  lateral  variations  in  the  crustal  and 
upper-mantle  P  wave  velocity  under  China. 

2.  Propagation  characteristics 
2.1  P  and  PP  waveforms 

Information  on  upper  mantle  structure  can  best  be  obtained  from  direct  P  phases  at  distances 
of  12.50  to  300.  In  the  12.50  to  20O  distance  range,  the  timing  and  amplitude  of  the  triplication 
arrivals  related  to  the  420  km  and  670  km  discontinuities  are  sensitive  to  velocity  gradients  in  the 
upper  mantle  (e.g.  Burdick,  1981).  At  epicentral  distances  between  20O  to  24o  these  triplication 
phases  arrive  in  a  very  short  time  window  and  are  not  very  diagnostic  of  velocity  variations  in  the 
upper  mantle.  Beyond  240  the  direct  P  phase  traverees  the  uppermost  mantle  at  a  steep  angle  and  is 
only  sensitive  to  the  average  lid  velocity  structure. 

Grand  &  Helmberger  (1984)  and  Lefevre  &  Helmberger  (1989)  have  demonstrated  that 
multibounce  P  and  S  waves  can  greatly  extend  the  distance  range  providing  information  on  the 
upper  mantle  velocity  structure.  Since  PP  triplication  phases  arrive  at  twice  the  distance  range  of 
corresponding  direct  P  wave  triplications  (25°  to  60°),  the  branches  associated  with  the  triplication 
arrivals  have  twice  the  time  separation  making  their  identification  much  easier.  In  addition, 
differences  between  P  and  PP  are  very  pronounced,  with  PP-P  differential  travel  times  and 
amplitudes  yielding  constraints  on  the  velocity  structure. 

Several  studies  have  concentrated  on  identifying  these  multibounce  triplication  phases  and 
examining  their  sensitivity  to  upper  mantle  velocity  structure  (Grand  &  Helmberger,  1985;  Lyon- 
Caen,  1986;  Lefevre  &  Helmberger,  1989;  Zhao  et  al.,  1991).  Using  synthetic  studies,  Schwartz 
and  Lay  (1993)  have  shown  that  the  presence  and  strength  of  precursors  to  the  PP  phase  is  highly 
diagnostic  of  velocity  gradients  in  the  lid.  For  positive  lid  velocity  gradients  strong  arrivals  are 
apparent  between  the  P  and  PP  wavetrains  at  distances  of  25°  to  34®,  with  less  energy  arriving  in 
this  time  window  for  constant  lid  velocities  and  only  weak  energy  for  negative  gradients.  These 
precursors  are  associated  with  energy  turning  in  the  lid  just  below  the  crust,  hugging  the  Moho 
(whispering  gallery  phases).  Beyond  34°  eticrgy  bottoming  deeper  in  the  upper  mantle  and 
multiply  reflecting  from  the  bottom  side  of  the  Moho  becomes  the  dominant  arrival  preceding  PP. 
This  arrival  is  stronger  for  negative  lid  velocity  gradients  which  force  energy  to  turn  deeper. 
Beyond  38°  little  precursory  energy  exists  other  than  interactions  with  the  transition  zone  (Neele  & 
Snieder,  1991).  At  distances  closer  than  25°  the  presence  of  precursors  to  PP  is  obscured  by  late  P 
triplication  arrivals  and  associated  coda.  The  whispering  gallery  phases  appear  as  impulsive 
arrivals  in  the  seismogram  because  much  of  their  low-frequency  energy  is  lost  through  tunneling  in 
the  lower  velocity  material  directly  beneath  the  lid  (Menke  &  Richards,  1980).  This  makes 
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whispering  gallery  phases  hard  to  discern  on  long-period  recordings.  The  presence  of  strong  lateral 
heterogeneity  may  also  obscure  the  identification  of  these  phases. 

Schwartz  &  Lay  ,  ■  993)  showed  that  crustal  reverberations  contribute  significantly  to  the  later 
part  of  the  PP  waveirain,  while  the  late  coda  of  PP  is  strongly  affected  by  surface  wave  scattering 
(Neele  &  Snieder,  1991).  PP-P  differential  travel  times  are  also  strongly  influenced  by  crustal 
thickness,  an  important  feature  that  can  easily  be  misinterpreted  as  variations  in  lid  velocity 
structure.  The  timing  of  the  PP  phase  can  also  be  confused  with  the  arrival  of  the  underside 
reflection  from  the  Moho.  Crustal  thicknesses  can,  to  some  extent,  be  interrogated  from  the 
characteristic  period  of  the  early  coda  trailing  the  main  P  and  PP  arrivals,  which  is  dominated  by 
crustal  reverberations.  Using  this  information  is  difficult  in  the  presence  of  strong  lateral  variations 
in  crustal  thickness,  since  for  the  P  and  PP  phases  these  reverberation  are  generated  at  two  or  three 
locations  along  the  earth's  surface,  respectively.  The  timing  between  the  underside  reflection  off 
the  Moho  and  the  first  energy  of  the  PP  phase  puts  another  constraint  on  the  crustal  thickness  at  the 
midpoint  of  the  raypath.  Because  the  PP  wave  has  distortion  due  to  being  a  maximum  time  phase, 
it  necessary  to  make  synthetic  seismograms  to  model  the  time  differences.  Since  there  are  so  many 
factors  contributing  to  the  timing  of  the  PP  phase  it  is  often  necessary  to  use  a  priori  velocity 
structure  information.  For  many  regions  detailed  information  on  crustal  thicknesses  is  available 
from  other  (e.g.  refraction)  studies.  Constraining  the  crustal  thicknesses  in  this  way  allows 
properties  of  the  upper  mande  to  be  more  reliably  determined. 

2.2  Regional  phases 

The  crustal  and  upper  mande  characteristics  that  we  attempt  to  model  using  the  P  and  PP 
waveforms  have  a  strong  influence  on  the  propagation  of  regional  phases.  The  Pn  and  Lg  phases 
in  particular  are  often  used  to  discriminate  between  earthquake  and  explosion  sources  and  to 
estimate  yields.  For  this  purpose  it  is  very  important  to  understand  how  the  excitation  and 
propagation  of  regional  phases  is  influenced  by  the  local  velocity  structure  in  order  to  enable 
correction  for  propagation  effects  along  the  raypath. 

The  Lg  phase  is  commonly  used  for  discrimination  purposes  because  of  its  stable  propagation 
properties  which  have  been  thoroughly  investigated  (e.g.  Hansen  et  al.,  1990;  Ringdal,  1991; 
Israelsson,  1992).  The  excitation  of  this  phase,  which  strongly  depends  on  how  much  S  wave 
energy  can  be  trapped  in  the  crust,  is  very  complicated.  For  an  explosion  source,  the  effectiveness 
of  Lg  excitation  is  governed  by  the  P  wave  velocity  near  the  source  and  the  S  wave  velocity  in  the 
uppermost  mantle  (Frankel,  1989).  If  the  P  wave  velocity  at  the  source  is  lower  than  the  Sn 
velocity  the  pS  phase  can  be  trapped  in  the  crust  by  (post-)critical  reflection  at  the  Moho.  In  this 
case,  the  pS  phase  is  the  dominant  energy  in  the  Lg  wavetrain.  The  amount  of  energy  that  can  be 
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trapped  in  the  crust  decreases  as  the  ratio  of  the  P  wave  velocity  near  the  source  and  the  S  wave 
velocity  in  the  uppermost  mantle  approaches  1.  This  can  be  attributed  to  the  fact  that  the  slowness 
range  (takeoff  angles)  that  will  lead  to  a  critical  reflection  at  the  Moho  narrows  as  this  ratio 
increases  (Xie  &  Lay,  1993).  A  similar  observation  can  be  made  for  sources  that  also  radiate  S 
wave  energy.  The  respective  criteria  for  P  and  S  energy  generated  near  the  source  to  be  trapped  in 
the  crust  are: 


sin  ip  > 
sin  is  > 


^Psource 

^Smantle 

^Ssource 

^Smantle 


(1) 
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where  ip  and  is  are  the  take-off  angles  of  P  and  S  waves.  If  the  P  wave  velocity  near  the  source  is 
greater  than  the  S  wave  velocity  in  the  uppermost  mantle  pS  cannot  be  trapped  in  the  crust  and 
energy  will  be  lost  through  radiation  to  the  mantle  at  each  consecutive  interaction  with  the  Moho.  In 
this  case  the  non-geometrical  S*  phase  will  be  the  main  contributor  to  the  Lg  phase  (Gutowski  et 
al.,  1984;  Frankel,  1989).  The  propagation  of  Lg  mainly  depends  on  the  thickness  and  structure  of 
the  crust;  a  thicker  crust  being  a  better  waveguide  for  this  phase  because  more  modes  can  be 
trapped. 

The  Pn  phase  is  critically  reflected  at  the  Moho  and  travels  along  the  bottomside  of  the  crust 
The  critical  reflection  of  Pn  energy  at  the  Moho  depends  on  the  P  wave  velocity  contrast  directly 
below  and  above  this  interface.  The  amount  of  energy  of  this  phase  that  can  leak  into  the  deeper 
mantle  depends  on  the  lid  velocity  gradients.  Higher  gradients  retain  more  P  energy  trapped 
directly  below  the  crust,  resulting  in  higher  Pn  amplitudes.  The  propagation  of  Pn  is  thus  strongly 
affected  by  tne  P  wave  velocity  structure  in  the  crust  and  uppermost  mantle. 


3.  Data  and  modeling  procedure 


In  order  to  obtain  P  and  PP  waveforms  suitable  for  modeling  of  the  crustal  and  upper  mantle 
structure  under  China  we  searched  the  Harvard-CMT  catalog  for  events  with  mb  larger  than  5.2  in 
the  southeast  Asian  region  which  were  recorded  within  a  14°  to  40®  epicentral  distance  range  by  at 
least  2  of  the  CDSN  stations.  Data  obtained  from  the  IRIS  database  in  Seattle  were  selected  on 
quality  of  the  P  waveforms  and  subsequently  deconvolved  by  the  instrument  response  to  yield 
displacement  traces.  We  obtained  recordings  from  6  broadband  seismic  stations  of  the  CDSN  array 
which  are  mainly  located  in  the  northern  half  of  China  (Figure  2).  Of  the  54  events  examined,  we 
selected  13  earthquakes  with  simple  waveforms  and  well  determined  focal  mechanisms.  This  large 
data  reduction  is  due  to  the  requirement  of  having  high  quality  data,  necessary  for  a  successful 
modeling  effort.  Without  clearly  identifiable  phases,  waveform  modeling  is  likely  to  fail  since 
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different  characteristics  of  the  data  are  easily  confused.  Most  epicenters  are  located  in  the 
tectonically  active  regions  of  western  and  southern  China.  There  is  a  total  of  45  paths,  primarily 
sampling  western,  central  and  northern  China  with  virtually  no  coverage  in  the  southeast  (Figure 
2).  For  the  post- 1988  period  of  our  event  search  no  high  quality  data  in  southeast  Asia  were  found 
to  eliminate  the  lack  of  coverage  in  this  region.  It  will  therefore  be  difficult  to  put  constraints  on  the 
velocity  structure  of  the  eastern  stable  platforms.  Our  models  are  primarily  for  the  active  fold  belts 
in  southwestern  China.  Table  1  gives  the  source  charactenstics  of  the  earthquakes  used  in  this 
study. 

We  match  the  observed  traces  by  one  dimensional  forward  modeling.  Our  approach  is  to  first 
model  traces  filtered  with  a  long-period  WWSSN  instrument  response.  This  filtering  facilitates 
determination  of  the  source  parameters  and  simplifies  identification  of  the  triplication  phases.  The 
objective  of  this  approach  is  to  obtain  a  first-order  velocity  model  that  matches,  to  the  extent 
possible,  all  long-period  waveforms.  This  model  then  can  serve  as  a  starting  model  for  matching 
the  broadband  data.  Modeling  of  the  broadband  data  refines  the  velocity  structure  information  and 
resolves  lateral  velocity  variations. 

Focal  mechanisms,  source  depths  and  source  time  functions  are  constrained  by  modeling  P 
waveforms  at  teleseismic  distances.  As  starting  parameters  we  use  the  Harvard  Centroid  Moment 
Tensor  (CMT)  solutions.  Our  events  range  in  magnitude  from  mb  =  5.5  to  6.8  with  source 
durations  of  up  to  7  seconds  (Table  1).  We  use  either  trapezoids  or  half  cycle  sinusoids  to  model 
the  source  time  functions.  For  events  with  magnitudes  higher  than  about  mb  =  7.0  such  a  simple 
representation  is  often  inadequate  and  the  more  complex  source  time  functions  obscure  the  velocity 
structure  information.  Events  smaller  than  about  mb  =  5.5  suffer  from  poorly  constrained  focal 
mechanisms  and  a  lower  signal  to  noise  ratio.  For  the  events  modeled,  source  complexity  and 
signal  to  noise  issues  mainly  complicate  matching  the  broadband  data;  long-period  waveforms  are 
easier  to  match. 

Differences  between  existing  procedures  for  generating  synthetic  seismograms  have  been 
thoroughly  investigated  and  are  well  understood  (e.g.  Burdick  &  Orcutt,  1979;  Chapman  &.  Orcutt, 
1985).  Because  the  WKBJ  method  of  Chapman  (1978)  is  exceptionally  fast  relative  to  other 
synthetic  seismogram  codes  it  is  well  suited  for  determining  the  source  parameters  from  teleseismic 
P  waves.  These  can  be  adequately  matched  with  the  inclusion  of  only  the  P,  pP  and  sP  ray  paths. 
Though  WKBJ  can  model  PP-P  differential  travel-times  and  gross  PP  amplifications,  it  fails  in 
matching  PP  waveform  characteristics  at  upper  mantle  distances  when  the  standard  practice  of 
including  only  primary  ray  paths  is  followed  (Schwartz  &  Lay,  1993).  Adding  enough  phases  to 
adequately  match  observed  PP  characteristics  is  a  time  consuming  and  slowly  converging  process. 
The  reflectivity  technique  developed  by  Fuchs  &  Muller  (1971)  completely  accounts  for  all  P-SV 
interactions  in  a  multi-layered  earth  and  can  adequately  match  observed  PP  waveforms.  After 
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determining  source  parameters  using  WKBJ  modeling  of  teleseismic  P  waveforms,  we  therefore 
used  reflectivity  synthetics  calculated  on  local  SUN-sparc2  stations  and  CRAY  or  Stellar 
computers  at  Lawrence  Livermore  National  Laboratory  (LLNL)  for  modeling  of  the  complete  P 
and  PP  waveforms  to  obtain  velocity  structure  information. 

4.  Results 

4.1  Long-period  modeling 

Our  first  modeling  objective  is  to  establish  a  homogeneous  reference  velocity  model  which 
can  match  observed  long-period  waveforms  on  average  for  the  whole  of  China.  In  order  to  assess 
what  range  of  upper  mantle  velocities  to  expect  for  China,  long-period  synthetics  were  initially 
generated  for  a  suite  of  3  velocity  models  (Figure  3)’ranging  from  the  fast  shield-type  models  S25 
for  the  Canadian  shield  (Lefevre  &  Helmberger,  1989)  and  K8  for  northwestern  Eurasia  (Given  & 
Helmberger,  1980),  to  the  slow  T7  model  (Burdick  &  Helmberger,  1978)  for  the  active  tectonic 
region  of  the  western  United  States.  Figure  4  shows  P  and  PP  waveforms  at  distances  ranging 
from  28^  to  36°  for  paths  through  different  tectonic  regions  of  China.  The  data  and  synthetic  traces 
are  aligned  on  their  first  motions,  thus  we  do  not  consider  absolute  travel  times.  This  makes  the 
matching  of  observed  waveforms  by  synthetics  less  sensitive  to  the  exact  event  locations,  thereby 
reducing  the  unknowns  that  have  to  be  considered  in  the  modeling  process.  Data  and  synthetics  are 
normalized  relative  to  their  peak  amplitudes.  Comparison  of  synthetics  generated  for  models  K8 
and  T7  with  the  data  shows  that  PP  arrives  too  early  for  K8  and  too  late  for  T7.  This  proves  to  be 
generally  true  for  observations  in  China.  At  distances  closer  than  about  30°  few  clear  PP  arrivals 
can  be  seen  in  the  long-period  data.  PP-P  differential  travel  times  are  strongly  influenced  by  P 
wave  velocities  in  the  uppermost  200  km  of  the  earth,  the  depth  range  where  differences  between 
models  K8  and  T7  are  most  pronounced. 

In  order  to  more  reliably  determine  P  wave  velocities  in  the  lid  we  constrain  the  crustal 
structure  by  using  results  of  previous  studies.  Several  refraction  studies  (Yuan  et  al.,  1986; 
Mooney,  1992)  have  shown  that  crustal  thickness  exceeds  45  km  for  most  of  China,  with 
maximum  thickness  of  up  to  70  km  for  the  Tibetan  plateau.  Since  most  paths  used  in  this  study 
cross  central  and  western  China,  a  50  km  thick  crust  for  our  preferred  model  WCH  (west  China)  is 
a  good  estimate  to  use  for  a  first  order  velocity  profile  suitable  for  this  region  (Figure  3,  Table  2). 
Holt  &  Wallace  (1990)  find  Pn  velocities  for  southern  China  ranging  from  8.0  km/s  for  the 
southeast  to  8.25  km/s  for  the  Tibetan  plateau.  Model  WCH  has  a  Pn  velocity  of  8.15  km/s  and  a 
low  velocity  zone  between  100  km  and  175  km,  necessary  to  slow  down  the  first  arrival  of  the  PP 
wave  train.  Figure  4  shows  that  model  WCH  is  able  to  explain  PP-P  differential  travel  times  over  a 
wide  range  of  distances  for  markedly  different  tectonic  regions. 
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The  synthetics  generated  for  model  K8  show  strong  PP  coda,  a  feature  that  is  not  observed  in 
the  data.  This  coda  mainly  consists  of  crustal  reverberations  (Schwartz  &  Lay,  1993),  the  strength 
of  which  is  determined  by  P  wave  velocities  directly  below  and  above  the  Moho.  The  reflection 
coefficient  at  the  Moho  increases  for  higher  Pn/Pcrust  ratios,  thus  increasing  the  energy  of  the 
crustal  bounces.  Since  we  adopted  the  same  crustal  velocity  of  6.4  km/s  as  in  model  K8  the  mantle 
lid  velocity  for  our  first  order  model  WCH  should  be  lower  than  that  of  model  K8.  The  Pn  velocity 
of  8.15  km/s  for  model  WCH  predicts  the  observed  PP  coda  well.  Also  apparent  from  Figure  4  is 
that  model  T7  predicts  too  much  energy  in  the  later  part  of  the  PP  arrival  itself.  The  amplitude  of 
the  PP  transition  zone  triplication  phases,  relative  to  shallow  turning  PP  energy,  is  strongly 
influenced  by  P  wave  velocities  in  the  lower  lid.  Lower  velocities  in  this  depth  range  tend  to 
decrease  the  amplitude  of  the  early  part  of  the  PP  wavetrain  for  which  the  energy  traverses  the 
lower  lid  at  a  shallow  angle.  This  indicates  that  our  first  order  model  needs  higher  velocities  than 
those  of  model  T7  in  the  lower  lid.  The  less  pronounced  low  velocity  zone  of  model  WCH  indeed 
better  predicts  relative  amplitudes  of  the  PP  triplication  phases  (Figure  4).  Though  difficult  to 
discern  in  the  long-period  data,  the  energy  arriving  in  the  window  between  P  and  PP  seems  well 
modeled;  another  indication  that  velocity  gradients  in  the  WCH  lid  are  suitable  for  the  China 
region.  Velocities  between  200  km  and  the  420  km  discontinuity  of  model  WCH  are  similar  to 
those  of  other  models  (Figure  3)  and  are  constrained  by  differential  travel  times  of  the  PP 
triplication  phases  and  overall  PP  amplitudes  (Figure  4).  The  effect  of  increasing  velocities  in  this 
depth  range  is  to  decrease  differential  travel  times  of  the  PP  triplication  phases  and  to  increase  PP/P 
amplitude  ratios.  These  features  are  in  general  well  predicted  by  model  WCH. 

A  P  and  PP  travel-time  curve  for  model  WCH  is  shown  in  Figure  5.  Figure  6  illustrates  that 
model  WCH  can  also  match  observed  P  triplication  behavior  for  paths  throughout  China. 
Comparison  with  Figure  4  shows  that  waveform  differences  at  P  triplication  distances  predicted  by 
the  three  models  are  considerably  smaller  than  PP  triplication  differences.  Beyond  18°,  differences 
in  the  synthetic  waveforms  for  the  different  models  are  very  small.  Given  the  sparse  data  coverage, 
adding  multibounce  phases  in  the  modeling  effort  is  thus  vital  for  a  detailed  study  of  the  upper 
mantle  velocity  structure.  Although  individual  features  of  the  data  may  be  better  matched  by  one  of 
the  slower  or  faster  models,  overall  WCH  best  explains  the  data.  It  is  remarkable  that  in  this  very 
heterogeneous  region  complex  long-period  waveform  characteristics  for  a  large  number  of  paths 
can,  to  a  great  extent,  be  explained  by  a  single  velocity  model. 

Although  model  WCH  can  match  observed  waveforms  for  paths  throughout  China  (Figures 
2,4  &  6)  it  is  evidently  more  characteristic  for  the  active  tectonic  region  of  southwestern  China  with 
its  thick  crust  and  fairly  high  Pn  velocity,  and  lateral  deviations  from  this  average  model  may  be 
expected.  This  first  order  model  has  a  number  of  remarkable  features  making  it  unusual  in 
character.  It  has  a  relatively  thick  crust  (50  km),  a  feature  that  contrasts  with  what  is  commonly 
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observed  for  active  tectonic  regions  (Mooney  &  Braille,  1989).  Model  WCH  also  has  a  Pn  velocity 
(8.15  km/s)  that  is  distinctly  higher  than  normally  seen  for  active  tectonic  regions.  The  lid  velocity 
structure  is  not  typical  for  an  active  tectonic  region,  or  for  a  stable  shield,  but  shows  an 
intermediate  behavior. 

Indications  of  the  range  of  variations  in  lid  velocity  structure  we  might  expect  beneath  China 
can  be  seen  from  the  long-period  waveforms  shown  in  Figure  7.  The  data  for  event  2,  recorded  at 
station  KMI  is  best  matched  by  model  T7.  The  considerably  slower  crust  and  lid  of  this  model 
pulls  apart  the  P  triplication  arrivals  enough  to  correctly  predict  the  arrival  of  energy  turning 
between  the  420  and  670  km  discontinuities  (CD  branch)  relative  to  energy  turning  above  the  420 
km  discontinuity  (AB  branch).  The  path  from  event  4  to  station  BJI  is  more  complicated.  The  CD 
branch  arrives  shortly  after  the  AB  branch,  a  feature  that  is  best  predicted  by  model  S25.  Energy 
turning  below  the  670  km  discontinuity  (EF  branch)  is  relatively  slow  which  is  best  matched  by 
model  WCH.  The  shape  and  amplitude  of  the  AB  branch  is  best  modeled  by  K8.  Since  the  timing 
of  energy  turning  between  the  420  and  670  km  discontinuities  is  influenced  more  by  the  shallow 
velocity  structure  than  energy  turning  below  the  670  km  discontinuity  this  indicates  that  velocities 
in  the  lid  should  be  higher  than  those  of  model  WCH  and  that  velocities  at  larger  depths  should  be 
somewhat  lower  than  those  of  model  WCH  to  compensate  for  the  effect  of  the  higher  lid  velocities 
on  the  EF  branch.  Similar  observations  can  be  made  for  the  raypath  from  event  4  to  station  WMQ. 
The  P  phase  for  this  path  is  best  explained  by  model  S25  whereas  the  arrival  of  PP  is  best 
predicted  by  model  T7,  although  the  PP  amplimde  is  overestimated  by  this  model.  The  best  match 
of  the  synthetic  for  model  K8  to  the  data  of  event  8,  observed  at  station  BJI  is  explained  by  this 
path  sampling  basin  type  material  over  its  whole  extent  (Figure  2).  In  summary,  of  the  45  paths  we 
used  in  the  long-period  modeling  process,  22  where  best  explained  by  model  WCH,  5  paths  were 
clearly  faster  and  1  path  was  slower,  the  remaining  paths  revealing  no  clear  preference  for  any  of 
the  models  considered.  An  overview  of  the  complete  long-period  modeling  results  is  given  in 
appendix  A.  Faster  and  slower  paths  are  not  confined  tc  distinct  regions,  but  are  intermixed, 
indicating  that  the  structure  of  the  crust  and  upper  mantle  under  China  cannot  simply  be  subdivided 
into  typical  shield  type  (S25,  K8)  regions  or  typical  active  tectonic  regions  (T7)  but  instead  shows 
a  more  complex  intermediate  behavior.  Using  the  broadband  waveforms  we  will  try  to  discriminate 
regionalized  deviations  from  the  average  model  WCH. 

4.2  Broadband  modeling 

Equipped  with  a  good  model  for  the  average  lid  velocity  we  can  use  precursors  to  PP  to 
interrogate  velocity  gradients  in  the  upper  mantle  directly  below  the  crust.  Energy  in  the  window 
between  P  and  PP  is  most  sensitive  to  these  gradients  and  will  therefore  suffer  least  from  non- 


10 


uniqueness  of  the  velocity  information  we  can  extract  from  the  data.  Two  examples  of  how  we  can 
use  precursors  to  PP  to  determine  the  lid  velocity  gradients  are  shown  in  Figure  8.  Model  WCH 
predicts  too  much  energy  preceding  PP  relative  to  the  observed  waveform  of  event  12  at  station 
BJI.  By  flattening  the  low  velocity  zone  this  energy  is  effectively  reduced.  The  reason  for  the 
reduced  velocities  in  the  200-400  km  depth  range  is  to  optimize  the  timing  of  the  PP  triplication 
arrivals,  the  thinner  crust  better  matches  the  timing  between  energy  reflected  off  the  Moho  and 
energy  reflected  off  the  crust.  The  WCH  synthetic  for  the  same  event  observed  at  station  HIA 
underpredicts  the  energy  arriving  directly  before  PP.  This  amplitude  car  be  increased  by  making 
the  low  velocity  zone  more  pronounced,  thus  increasing  the  velocity  gradients  at  both  the  top  and 
bottom  of  the  low  velocity  zone.  Considerably  higher  gradients  are  needed  to  significantly  increase 
the  precursory  energy.  Since  model  WCH  correctly  predicted  the  timing  of  the  PP  arrival,  the 
average  velocity  in  the  lid  needs  to  remain  the  same,  requiring  an  increased  Pn  velocity  Thus  by 
changing  the  shape  of  the  low  velocity  zone  we  can  control  the  amount  of  energy  arriving  in  the 
time  window  between  the  P  and  PP  phases.  Increasing  the  negative  velocity  gradient  causes  mote 
energy  to  turn  deeper  in  the  manfle  in  the  positive  velocity  gradient  below  the  low  velocity  zone. 
This  leads  us  to  conclude  that,  even  at  distances  as  close  as  31°,  energy  turning  below  the  lid  and 
multiply  reflecting  off  the  bottom  side  of  the  Moho  is  the  dominant  arrival  in  this  time  window. 

More  evidence  for  deviations  from  model  WCH  in  the  southwestern  part  of  China  can  be 
seen  for  other  recordings  of  event  12  (Figure  9).  At  about  15°  the  P  triplication  arrivals  are  clearly 
separated  (Figure  5),  with  the  AB  branch  arriving  first  and  the  EF  branch  arriving  last.  Near  22° 
the  AB  branch  almost  coincides  with  the  CD  branch  for  model  WCH,  both  arriving  about  7 
seconds  before  EF.  The  AB  and  CD  branches  together  form  the  first  upswing  in  the  LZH 
seismogram  followed  by  the  CD  depth  phases  which  can  be  seen  as  distinct  negative  arrivals,  the 
EF  branch  being  only  a  weak  arrival  in  the  seismogram.  The  P  triplication  phases  as  predicted  by 
Hi  lei  WCH  arrive  in  a  too  short  time  window  for  both  waveforms  shown  in  Figure  9.  Lowering 
the  velocities  at  depths  greater  than  200  km  effectively  pulls  apart  the  triplication  arrivals.  A 
secondary  effect  of  these  lower  velocities  on  the  synthetic  waveform  for  station  WMQ  is  that  the 
energy  of  the  AB  branch  is  reduced.  This  can  be  counteracted  by  increasing  velocities  in  the  lid. 
These  higher  lid  velocities  only  slightly  affect  the  timing  of  the  triplication  arrivals  because  the  P 
energy  traverses  this  part  of  the  mantle  at  a  steep  angle.  The  added  crustal  thickness  for  this  path 
gives  a  good  fit  to  the  observed  P  wave  coda.  The  increased  Pn  velocity  for  the  waveform  recorded 
at  station  LZH  optimizes  the  timing  between  energy  turning  above  the  420  km  discontinuity  and  the 
depth  phases  associated  with  energy  turning  between  the  420  and  670  km  discontinuities. 

The  velocity  structure  of  the  stable  eastern  platforms  can  be  interrogated  to  a  limited  extent 
using  paths  branching  to  the  northeast  from  event  3,  partially  traversing  these  platforms  (Figure  2). 
The  strongest  indications  for  the  lid  velocity  structure  can  again  be  derived  from  modeling  PP 
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characteristics.  For  both  waveforms  shown  in  Figure  10  the  timing  and  general  shape  of  the  PP 
phase  can  best  be  explained  by  a  constant  velocity  lid  and  crustal  thicknesses  that  are  considerably 
smaller  than  50  km.  For  the  path  to  station  MDJ  a  somewhat  lower  Pn  velocity  (8.12  km/s)  is 
needed  to  match  the  timing  of  PP  energy.  The  thin  crust  (27  km)  for  this  path  also  matches  the 
crustal  reverberations  trailing  the  PP  wavetrain. 

Event  4  was  used  to  investigate  the  crustal  and  upper  mantle  velocity  strucmre  in  the  fold  belt 
region  of  northern  China  and  Mongolia.  The  waveforms  shown  in  Figure  1 1  reveal  evidence  for 
lower  Pn  velocities  (8.0  km/s)  combined  with  a  positive  velocity  gradient  in  the  lid.  For  the  path 
from  event  4  to  station  BJI  model  WCH  predicts  the  second  (CD)  branch  of  the  triplication  arrivals 
to  arrive  too  close  after  the  first  (AB)  branch.  The  effect  of  the  lower  velocities  in  the  200-420  km 
depth  interval  is  to  pull  apart  these  triplication  arrivals.  Higher  velocities  in  the  lower  lid  increase 
the  amplitude  of  the  CD  relative  to  the  AB  branch.  However,  these  higher  velocities  also  shorten 
the  time  window  of  the  triplication  arrivals.  The  average  lid  velocity  needs  to  remain  equal  to  that 
of  model  WCH,  requiring  a  lower  Pn  velocity  and  thus  a  positive  gradient  in  the  lid.  The  effect  of 
this  gradient  on  the  synthetic  waveform  for  station  WMQ  is  to  better  predict  the  shape  of  the  first 
positive  upswing  in  the  P  arrival.  A  lid  with  a  positive  velocity  gradient  also  better  matches  the 
observed  arrival  of  the  second  (CD)  branch  of  the  PP  triplication  which  was  predicted  too  late  for 
model  WCH.  The  second  negative  peak  of  this  branch,  a  depth  phase,  is  still  poorly  predicted 
however.  We  did  not  try  to  improve  the  match  because  it  might  also  be  caused  by  a  slightly 
erroneous  focal  mechanism.  This  is  not  unlikely  since  the  amplitude  overprediction  of  the  onset  of 
the  P  phase  can  also  be  attributed  to  a  focal  mechanism  error.  We  infer  this  because  the  mismatch  is 
consistently  present  for  synthetics  generated  with  velocity  models  as  different  as  S25  and  T7.  The 
onset  of  PP  arriving  shortly  after  P  and  energy  in  between  these  two  phases  does  seem  better 
matched  by  the  gradient  lid  model. 

We  have  found  significant  deviations  in  crustal  and  upper  mantle  structure  from  the  average 
China  model  WCH.  Figure  12  shows  our  preferred  models  for  the  Tibetan  plateau  (TP),  the 
eastern  platforms  (CHE)  and  northern  China  (NCH)  together  with  the  average  model  WCH  and  P 
and  S  wave  models  of  previous  studies;  Figure  13  summarizes  lateral  variations  in  the  upper 
mantle  velocity  structure  under  China  at  different  depth  intervals.  While  the  model  variations  are 
not  unique  we  emphasize  the  more  robust  features.  An  overview  of  the  modeling  results  is  given  in 
appendix  B.  A  total  of  6  events  (19  paths)  were  used  for  broadband  modeling.  Source  mechanisms 
of  the  other  events  could  not  be  constrained  well  enough  to  confidently  natch  the  broadband  data. 
In  contrast  to  the  long-period  synthetics  which  are  not  very  sensitive  to  small  variations  in  source 
mechanism,  details  of  the  broadband  ^ j  .tic  waveforms  can  change  significantly. 

Though  absolute  lid  velocities  of  the  Russian  P  wave  model  CE200  (Figure  12)  of  Gamero  et 
al.  (1992)  are  higher  than  those  of  model  WCH,  the  lid  velocity  gradients  of  the  two  models  are 
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very  similar.  The  most  striking  difference  is  the  220  km  discontinuity  their  model  indicates.  We 
included  a  220  km  discontinuity  in  our  models  to  see  whether  the  presence  of  such  a  velocity  jump 
was  required  by  our  data.  However,  a  velocity  contrast  similar  to  that  of  model  CE200  changed 
only  minor  details  of  the  broadband  synthetics.  We  therefore  cannot  discern  whether  or  not  a  220 
km  discontinuity  is  present  in  the  upper  mantle  under  China.  The  presence  of  multilayering  and 
lateral  heterogeneity  in  the  crust  has  been  well  documented  (e.g.  Mangino  &  Ebel,  1992;  Mooney, 

1992) .  We  investigated  the  effect  of  a  multilayered  crust  on  the  broadband  synthetics  and  found 
that  we  could  not  confidently  constrain  depth  dependent  velocity  information  in  the  crust  with  our 
dataset.  Furthermore,  no  evidence  for  lateral  variations  in  crustal  velocity  was  observed,  although 
S  wave  velocities  in  Tibet  are  believed  to  be  anomalously  low  (F.T.  Wu,  personal  communication, 

1993) . 

Figure  13a  shows  that  crustal  thicknesses  average  50  km  in  the  Tibetan  region  with  2  paths 
indicating  thicknesses  exceeding  this  value.  Somewhat  lower  crustal  thicknesses  (43-50  km)  are 
found  for  the  northern  fold  system  and  considerably  lower  thicknesses  (27-35  km)  for  3  paths 
partly  traversing  the  eastern  platforms.  These  values  are  somewhat  lower  than  results  of  previous 
studies  (e.g.  Yuan  et  al.,  1986;  Mooney,  1992),  but  we  have  limited  resolution  of  crustal 
thickness.  A  study  of  the  receiver  structure  beneath  the  6  CDSN  stations  by  Mangino  &  Ebel 
(1992)  revealed  10-20  km  thick  transition  zones  between  the  crust  and  upper  mantle,  with 
velocities  sharply  increasing  from  about  6.0-6.4  km/s  to  a  fixed  preassumed  local  Pn  velocity,  for 
a  number  of  these  stations.  Such  a  transition  zone  cannot  be  resolved  by  our  data  and  its  existence 
might  be  responsible  for  the  shallower  Moho  depths  of  this  study.  Average  crustal  velocities  of  6.4 
km/s  are  suitable  for  all  the  paths  used  in  this  study. 

We  find  Pn  velocities  (Figure  13b)  ranging  from  8.15-8.25  km/s  for  the  Tibetan  plateau  and 
as  low  as  8.0  km/s  for  the  northern  fold  belt.  Paths  traversing  northeastern  China  need  Pn 
velocities  of  about  8.12-8. 14  km/s.  Our  Pn  velocity  range  found  for  the  Tibetan  plateau  agrees  well 
with  results  of  previous  studies  (e.g.  Holt  &  Wallace,  1990;  Mooney,  1992).  Highest  Pn  velocities 
(8.20-8.25  km/s)  can  be  observed  for  path  originating  in  the  western  part  of  the  plateau, 
compatible  with  the  unusually  high  Rayleigh  wave  phase  velocities  found  by  Brandon  & 
Romanowicz  (1986)  for  this  part  of  Tibet.  Refraction  studies  indicate  Pn  velocities  averaging  8.0 
km/s  near  BJI  and  values  up  to  8.2  km/s  north  of  KMI,  a  velocity  range  that  is  in  agreement  with 
Pn  velocities  we  find  for  the  easternmost  paths.  Besides  a  single  refraction  profile  in  the  Tarim 
basin  which  indicates  a  Pn  velocity  of  7.9  km/s  (Mooney,  1992)  no  values  have  been  reported  for 
the  northern  part  of  China  and  Mongolia. 

This  study  indicates  that  the  Tibetan  plateau  has  a  thin  lid  underlain  by  a  low  velocity  zone 
with  lower  velocities  extending  below  200  km  depth  (Figure  13c,d).  This  observation  is  very 
similar  to  S  wave  velocity  results  (Figure  12)  of  Lyon-Caen  (1986)  and  Wang  et  al.  (1989,  not 
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shown).  The  low  velocity  zone  in  our  model  TP  is  located  somewhat  deeper  than  in  Lyon-Caen's 
(1986)  model,  and  indicates  a  lid  of  about  50  km  thickness.  Given  previous  evidence  for  high  Pn 
velocities  under  the  Tibetan  plateau,  the  timing  of  PP  relative  to  P  for  observed  waveforms  of  paths 
traversing  this  region  requires  the  presence  of  a  distinct  low  velocity  zone  in  the  upper  mantle.  The 
sensitivity  of  precursors  to  PP  to  velocity  gradients  in  the  uppermost  mantle  puts  constraints  on  the 
Tibetan  lid  velocity  structure  and  thereby  also  on  the  location  of  the  low  velocity  zone  (Figure  8). 

Using  Pn/Pl  ratios  Holt  &  Wallace  (1990)  determined  positive  lid  velocity  gradients  of  0.18- 
0.25  km/s  per  100  km  with  a  100  km  thick  lid  for  southern  Tibet.  Although  our  sources  are  located 
in  their  region  of  study  their  result  is  not  necessarily  in  disagreement  with  what  we  find  since  most 
of  our  paths  extend  well  north  of  this  region.  Two  paths  from  our  southernmost  event  2  mainly 
traversing  the  Tibetan  plateau  indicate  either  no  (station  LZH)  or  a  shallow  low  velocity  zone 
(WMQ)  indicating  that  a  transition  zone  from  a  shield  type  lid  in  southernmost  Tibet,  as  inferred 
from  Holt  &  Wallace's  data,  to  a  weaker  lid  in  central  and  northern  Tibet  may  exist.  However,  due 
to  a  lack  of  paths  purely  traversing  southern  Tibet  we  cannot  clearly  resolve  this. 

From  our  results  we  infer  that  the  greater  part  of  the  Tibetan  plateau  is  presently  not  underlain 
by  a  shield-type  structure.  Thus,  if  the  whole  of  Tibet  were  originally  underthrust  by  the  Indian 
plate  this  observation  requires  that  the  cold  lithospheric  material  has  disappeared  in  a  relatively 
short  time  period.  The  high  Pn  and  Sn  velocities  observed  for  the  Tibetan  plateau  do  not 
necessarily  imply  the  underthrusting  of  cold  lithospheric  material.  Crustal  shortening  and 
thickening  requires  that  the  underlying  mantle  material  also  undergoes  shortening  and  thickening, 
moving  cold  isotherms  downward.  In  this  way  higher  velocities  in  the  uppermost  mantle  can  be 
produced  by  increased  pressures  due  to  the  added  crustal  thickness  in  the  absence  of  higher 
temperatures  (Black  &  Braille,  1982;  Holt  &  Wallace,  1990).  We  prefer  the  crustal  shortening 
model  for  the  evolution  of  the  Tibetan  plateau  because  it  is  able  to  explain  a  wide  range  of 
observations  on  the  crustal  and  upper  mantle  structure  in  this  region.  Underthrusting  seems  to  have 
only  played  a  role  in  the  shaping  of  the  southernmost  part  of  the  plateau  and  the  Himalayas 
(Molnar,  1988). 

Lowered  temperatures  in  the  lid  below  Tibet  will  lead  to  horizontal  temperature  gradients  in 
the  mantle.  These  are  likely  to  introduce  thermal  convection  with  horizontal  motions  directed 
toward  the  colder  material  beneath  central  Tibet  This  requires  upwelling  warmer  mantle  material  in 
the  neighboring  regions.  A  cartoon  illustrating  this  scenario  is  given  in  Figure  22  of  Molnar 
(1988).  A  recent  tomographic  study  of  lateral  variations  in  P  wave  velocity  beneath  the  Tibetan 
plateau  by  Zhao  &  Xie  (1993)  reveals  a  velocity  image  with  highs  and  lows  coinciding  with 
proposed  locations  of  upwelling  and  downwelling  in  the  evolution  model  of  Molnar  (1988). 
Surface  evidence  for  mantle  upwelling  can  be  seen  from  C^iatemary  volcanism  in  the  Chang  Thang 
region  in  north-central  Tibet.  A  shallower  crust  together  with  lower  Q  and  Sn  values  for  the  lid. 
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related  to  higher  upper  mantle  temperatures,  has  been  documented  (e.g.  Ni  &  Barazangi,  1983; 
Bourjot  &  Romanowicz,  1992)  for  this  region.  The  S  wave  velocity  model  of  Zhao  et  al.  (1991) 
with  its  higher  crustal  velocity  and  slower  lid  than  Lyon-Caen's  (1986)  in  which  a  low  velocity 
zone  is  virtually  absent  (Figure  12)  also  seems  appropriate  for  the  Chang  Thang  region.  While  we 
find  evidence  for  the  presence  of  lateral  heterogeneity  in  the  Tibetan  lid  (Figure  13a,b),  the  length 
and  orientation  of  our  paths  is  such  that  we  cannot  constrain  the  location  and  extent  of  such 
inhomogeneities.  Combining  lid  thicknesses  inferred  by  Lyon-Caen  (1986)  and  Holt  &  Wallace 
(1990)  with  ours,  indicates  a  northward  thinning  of  the  Tibetan  lid.  This  suggests  southward 
motion  of  warmer  mantle  material  from  the  Chang  Thang  region,  gradually  weakening  the  lid  from 
below.  Molnar  (1988)  theorizes  that  the  extremely  high  phase  velocities  found  for  westernmost 
Tibet  (Brandon  &  Romanowicz,  1986),  which  are  confirmed  by  high  Pn  velocities  we  find  for  that 
region,  might  indicate  active  crustal  shortening  accompanied  by  the  sinking  of  cold  material.  Our 
modeling  results  show  lower  than  average  velocities  between  200-420  km  depth  for  the  whole  of 
Tibet  (Figure  13d).  These  indicate  the  presence  of  warmer  mantle  material  in  a  broad  region  in  this 
depth  range.  Thus,  very  strong  lateral  heterogeneity  in  upper  mantle  temperatures  seems  to  exist 
below  the  Tibetan  plateau. 

Holt  &  Wallace  (1990)  find  strong  evidence  for  a  zero  gradient  model  with  a  100  km  thick  lid 
below  southeastern  China,  similar  to  what  our  models  infer  for  eastern  China  (model  CHE,  Figure 
12).  In  contrast,  Zhao  et  al.  (1991)  indicate  an  S  wave  model  with  a  pronounced  low  velocity  zone 
for  the  Yangtze  and  Sino-Korean  platforms.  However,  since  their  sources  are  located  at  the  edge  of 
the  Tibetan  plateau  and  their  stations  are  mostly  underlain  by  oceanic  lithosphere,  their  paths  cross 
three  very  different  tectonic  provinces.  A  P  wave  velocity  lid  similar  to  the  model  of  Zhao  et  al. 
(1991)  would  slow  down  the  PP  phase  too  much,  unless  such  a  lid  was  accompanied  by  a  higher 
Pn  velocity.  Our  Pn  velocity  range  of  8.12-8.14  km/s  for  the  easternmost  paths  is  in  agreement 
with  the  range  found  by  deep  refraction  studies  (Mooney,  1992).  Furthermore,  the  constant 
velocity  lid  predicts  hard  to  match  details  of  observed  broadband  PP  behavior  well  (Figure  10).  We 
therefore  infer  that  the  eastern  Chinese  platforms  can  be  characterized  by  a  thick  constant  velocity 
lid  indicating  a  stable  evolution  for  the  upper  mantle  in  this  region  in  recent  times. 

The  positive  lid  velocity  gradient  model  this  study  finds  for  northern  China  and  Mongolia  is 
preliminary  and  should  be  viewed  with  cate.  We  only  modeled  3  paths  traversing  this  region  and, 
given  the  uncertainty  in  the  mechanism  of  event  4,  the  evidence  for  a  positive  lid  velocity  gradient 
model  is  not  very  strong.  Further,  no  evidence  either  confirming  or  contradicting  this  observation 
has  been  found. 

4.3  Propagation  of  regional  phases 
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Having  established  regionalized  velocity  models  for  China  we  perform  a  synthetic  study  to 
assess  how  these  velocity  models  affect  the  propagation  of  regional  phases.  Pn  and  Lg  phases  are 
often  used  in  regional  discrimination  of  explosion  and  earthquake  sources.  We  will  therefore 
analyze  these  phases  using  synthetics  generated  for  the  regionalized  velocity  models  and  explore 
differences  in  propagation  for  the  different  models. 

The  reflectivity  technique  of  Fuchs  &  Muller  (1971)  was  used  to  generate  synthetic 
seismograms  at  regional  distances.  We  used  an  explosive  source  located  at  500  m  depth  and  an 
earthquake  source  at  5  km  depth,  both  with  a  source  duration  of  1.5  seconds.  We  calculated 
seismograms  at  distances  of  1°  to  10°  at  1°  intervals.  The  station  azimuth  for  the  earthquake 
source  is  taken  to  be  22.50.  Due  to  computational  limits  we  were  required  to  limit  the  maximum 
frequency  of  the  synthetic  seismograms  to  2  Hz.  Observed  regional  waveforms  are  known  to  have 
a  frequency  content  of  up  to  5  Hz.  Though  this  somewhat  limits  us  in  the  extent  to  which  we  can 
quantify  the  effect  of  the  regionalized  models  on  the  propagation  of  the  Pn  and  Lg  phases,  most  of 
their  spectral  energy  is  usually  concentrated  well  below  2  Hz.  In  order  to  effectively  generate  the 
Lg  phase  for  the  explosion  source  we  include  a  thin  crustal  layer  at  the  free  surface  with  a  low 
enough  velocity  to  capture  the  pS  phase  in  the  crust  (Frankel,  1989).  We  use  a  3  km  thick  layer 
with  a  P  wave  velocity  of  4.4  km/s.  The  presence  of  such  low  velocities  in  the  uppermost  crust  is 
well  documented  (e.g.  Mangino  &  Ebel,  1992)  but  in  general  cannot  be  resolved  by  P  and  PP 
waveform  modeling.  The  low  velocity  layer  does  not  have  a  significant  effect  on  the  P  and  PP 
waveforms  we  used  in  this  study  because  of  its  limited  thickness. 

We  generated  synthetics  for  model  WCH  and  models  typical  for  the  Tibetan  plateau  (TP)  and 
eastern  (CHE)  and  northern  (NCH)  China  (Figure  12).  Pn/Sn  velocity  ratios  in  the  lid  are  1.731 
for  all  models.  Synthetic  broadband  vertical  component  displacement  waveforms  for  model  WCH 
and  the  explosive  source  are  shown  in  Figure  14.  The  synthetics  exclude  phase  velocities  smaller 
than  3.3  km/s  to  eliminate  the  high  amplitude  Rayleigh  wave  arrivals  from  the  seismograms.  We 
define  several  windows  for  making  measurements  of  the  regional  phase  energy.  The  Pn  window  is 
a  constant  time  window  starting  just  ahead  of  the  onset  of  Pn  (we  let  the  start  time  depend  on  the 
velocity  model  used)  with  a  time  length  of  10  s;  the  Lg  phase  arrives  in  the  3.3-3.7  km/s  velocity 
range.  As  can  be  seen  from  Figure  14,  the  Pn  window  is  not  very  useful  at  distances  closer  than  4° 
because  other,  higher  energy  phases  (e.g.  Pg)  enter  the  window.  The  first  distinct  arrival  outside 
this  window  is  PnPn-  Pg  marks  the  onset  of  high  frequency  energy  in  between  the  Pn  and  Lg 
windows.  The  long-period  component  in  the  Lg  window  consists  of  low  order  Rayleigh  wave 
modes.  The  onset  of  Rayleigh  wave  energy  can  be  identified  at  a  phase  velocity  of  about  4.2  km/s. 
The  higher  frequency  components  superimposed  on  the  long-period  motions  are  the  actual  Lg 
energy. 
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The  long-period  Rayleigh  wave  signals  are  not  recorded  by  the  short-period  instruments 
which  are  most  commonly  used  to  record  regional  waveforms.  In  order  to  simulate  the  response  of 
such  an  instrument  we  convolve  the  synthetic  displacement  traces  with  the  response  of  the  Obninsk 
(OBN)  Kimos  instrument,  typical  of  instruments  operating  in  the  C.I.S..  The  effect  of  this  filtering 
in  the  frequency  domain  is  shown  in  Figure  15.  The  frequency  band  of  the  instrument  extends  to 
about  8  Hz,  leaving  the  high  frequency  component  of  the  synthetics  unchanged.  The  filtered 
waveforms  are  presented  in  Figure  16.  The  high  frequency  arrivals  between  the  Pn  and  Lg 
windows  have  now  become  the  largest  arrivals  in  the  seismograms.  The  long-period  Rayleigh 
wave  energy  has  disappeared  from  the  Lg  window  and  the  Pn  phase  has  become  more  impulsive 
and  complex.  Because  many  regional  waveforms  are  now  being  recorded  by  broadband 
instruments  we  investigate  propagation  characteristics  of  the  Pn  and  Lg  phases  for  both  the 
broadband  and  the  short-period  signals.  We  do  this  by  calculating  root  mean  square  amplitudes  for 
the  2  time  windows  as  a  function  of  epicentral  distance. 

Figure  17  shows  the  root  mean  square  amplitudes  in  the  Pn  and  Lg  windows  for  the  short- 
period  waveforms  generated  with  the  4  velocity  models  derived  in  this  study.  Amplitude 
differences  for  the  average  model  WCH  and  model  TP  are  very  small  for  both  Pn  and  Lg  indicating 
that  the  higher  Pn  velocity  of  the  latter  model  does  not  change  the  propagation  characteristics  of 
these  phases  significantly.  At  distances  closer  than  4®  the  Pn  window  is  contaminated  by  the 
inclusion  of  other  phases  in  the  rms  window  and  we  will  therefore  only  comment  on  propagation 
characteristics  at  greater  distances.  Between  4^  and  10°  the  increased  amplitude  of  the  Pn  phase  for 
model  CHE  relative  to  the  average  model  WCH  can  be  attributed  to  the  reduced  crustal  thickness. 
The  effect  of  the  lower  P  wave  velocity  in  the  uppermost  mantle  for  this  model  is  to  decrease  the 
Pn  amplitude  because  the  increase  in  ratio  between  the  source  P  wave  velocity  and  the  Pn  velocity 
narrows  the  range  of  takeoff  angles  that  will  give  a  critical  reflection  at  the  Moho  and  thus  more 
energy  will  be  lost  through  radiation  into  the  mantle.  Thus,  the  effect  of  changes  in  crustal 
thickness  appear  more  important  than  baseline  shifts  in  lid  velocity. 

We  would  expect  Pn  amplitudes  for  model  NCH  to  be  lower  than  for  model  CHE  because 
NCH  has  both  a  thicker  crust  and  a  lower  Pn  velocity.  However,  in  contrast  we  find  higher  Pn 
amplitudes  increasing  with  distance.  These  observations  are  explained  by  the  positive  velocity 
gradient  in  the  lid  for  model  NCH.  The  positive  lid  gradient  causes  less  Pn  energy  traveling  along 
the  bottom  side  of  the  Moho  to  be  lost  through  radiation  into  the  upper  mantle.  Another  effect  is 
that  energy  that  is  almost  critically  reflected  at  the  Moho  and  penetrates  into  the  lid  at  near 
horizontal  angles,  turns  at  very  shallow  depth  and  can  contribute  to  energy  in  the  Pn  window.  This 
latter  effect  increases  with  distance  accounting  for  the  gain  in  Pn  energy  with  distance.  This 
increase  seems  to  flatten  off  near  10^  indicating  that  effects  of  damping  and  geometrical  spreading 
are  becoming  equally  important.  Though  baseline  shifts  in  P  velocities  in  the  uppermost  mantle  do 
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not  influence  Pn  propagation  characteristics  very  much,  variations  in  lid  velocity  gradients  can 
significantly  change  observed  Pn  amplitudes. 

The  small  difference  in  Lg  energy  between  models  WCH  and  TP  and  the  significant  decrease 
for  models  NCH  and  CHE  indicates  that,  like  i or  Pn,  crustal  thickness  is  the  most  important  factor 
controlling  Lg  energy,  with  a  thicker  crust  being  a  better  waveguide  for  this  phase.  Variations  in  lid 
velocities  seem  to  have  only  a  minor  effect  on  the  propagation  of  Lg.  Although  some  scatter  is 
visible,  the  propagation  of  Lg  seems  very  stable  and  the  amplitude  decay  with  distance  can  be  well 
matched  by  a  simple  power  law  equation.  The  powers  of  the  distance  decay  rates  are  given  in  table 
3.  For  the  different  models  these  range  from  2.35  to  2.46.  This  range  is  in  agreement  with 
observations  by  Blandford  (1981)  who  finds  empirical  Lg  decay  rates  of  for  the  active  tectonic 
region  of  the  western  United  States  to  for  the  stable  eastern  United  States.  Figure  17  further 
shows  that  differences  in  propagation  for  the  different  models  can  change  Pn/Lg  ratios  by  as  much 
as  an  order  of  magnitude.  At  distances  larger  than  4°  Pn/Lg  ratios  increase  with  distance  because 
the  decay  of  Lg  with  distance  is  faster  than  that  of  Pn,  with  the  slope  of  the  increase  depending  on 
the  difference  in  the  decay  rates  of  the  2  phases. 

Figure  18  shows  identical  measurements  for  the  broadband  waveforms.  The  Pn  amplitude 
curves  for  this  case  are  much  simpler.  Except  for  model  NCH  they  can  all  be  characterized  by 
simple  power  law  decays  (Table  4).  Also,  differences  in  rms  Pn  amplitudes  for  the  different 
models  are  much  smaller  than  for  the  short-period  waveforms.  The  amplitude  decay  in  the  Lg 
window  can  again  be  characterized  by  power  law  equations  (Table  3)  with  decay  rates  being 
significantly  smaller  than  for  the  filtered  waveforms.  The  long-period  Rayleigh  wave  energy  that 
now  dominates  in  this  window  (Figure  14)  damps  out  less  rapidly  than  the  higher  frequency  Lg 
modes.  Pn/Lg  ratios  now  uniformly  decrease  with  distance,  except  for  model  NCH.  The  behavior 
of  this  decrease  is  fairly  simple  and  the  slope  differences  are  governed  by  variations  in  the  Lg 
decay  rates.  Depending  on  the  structure  of  the  crust  and  upper  mantle  Pn/Lg  ratios  for  the 
broadband  synthetics  vary  by  a  factor  5. 

The  short-period  waveform  behavior  for  the  earthquake  source  is  very  similar  to  that  of  the 
explosion  source  (Figures  17  «&  19),  overall  amplitudes  of  both  the  Pn  and  Lg  phase  are  down  by 
about  a  decade  but  relative  differences  between  the  different  models  are  quite  small.  The  most 
marked  difference  is  the  change  in  slope  of  the  Lg  decay  for  the  north  and  east  China  models 
(Table  3).  In  these  cases  the  P  wave  velocity  near  the  source  is  larger  than  the  S  wave  velocity  in 
the  uppermost  mantle  and  the  pS  phase  cannot  be  captured  in  the  crust.  Therefore  the  main 
contribution  to  the  Lg  phase  will  come  from  direct  excitation  of  S  wave  energy.  Because  the  Lg 
variations  between  models  WCH  and  TP  are  small  the  effect  of  velocity  changes  in  the  uppermost 
mantle  appear  to  not  be  very  important.  Crustal  thickness  is  the  major  factor  governing  the  strength 
of  Lg.  The  markedly  slower  Lg  amplitude  decay  for  model  NCH  can  be  attributed  to  S  wave 
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energy  that  is  not  critically  reflected  at  the  Moho  turning  in  the  uppermost  mantle  and  contributing 
to  the  energy  in  this  window.  Pn/Lg  ratios  beyond  4®  are  overall  somewhat  lower  for  the 
earthquake  source  than  those  for  the  explosion  source,  with  a  stronger  effect  for  the  east  and  north 
China  models. 

Comparison  of  Figures  18  and  20  shows  that  the  broadband  Pn  behavior  for  the  earthquake 
source  is  very  similar  to  that  of  the  explosion  source.  Overall,  broadband  Pn  depends  less  on 
velocity  structure  than  Lg.  The  low  order  Rayleigh  waves  which  dominate  the  broadband  Lg 
window  exhibit  a  complex  behavior.  At  close  distances  the  energy  in  this  window  is  lowest  for 
models  NCH  and  CHE  but  because  the  amplitude  decay  for  these  models  is  slower  than  that  for 
models  WCH  and  CHE  they  give  higher  amplitudes  at  distances  greater  than  7°.  Model 
dependence  of  Pn/Lg  ratios  in  this  case  is  very  small. 

For  discrimination  between  earthquakes  and  underground  nuclear  explosions  we  would  like 
the  Pn/Lg  variations  between  earthquake  and  explosion  sources  to  be  pronounced  and  to  have  a 
simple  predictable  behavior.  For  discrimination  purposes  Pn/Lg  ratios  are  often  expected  to  be 
higher  for  nuclear  explosions  than  for  earthquakes,  mainly  due  to  enriched  high  frequency  Pn 
energy  for  explosions  and  enriched  Lg  for  earthquakes.  However,  the  real  earth  is  not  that  simple. 
Indeed,  several  observational  studies  have  either  confirmed  higher  Pn/Lg  ratios  for  earthquake 
sources  (e.g.  Blandford,  1981,  Chan  et  al.,  1990)  or  found  that  Pn/Lg  ratios  have  much  more 
complicated  behavior  (e.g.  Lynnes  &  Baumstark,  1991).  Shown  in  Figure  21  is  the  quotient  of 
these  ratios  for  the  explosion  and  earthquake  synthetics.  Because  we  use  identical  paths  for  the 
earthquake  and  explosion  source  this  quotient  is  mainly  sensitive  to  source  radiation  differences. 
For  the  short-period  frequency  band  the  explosion/earthquake  quotient  ranges  roughly  between  1 
and  3  at  distances  larger  than  40.  While  we  find  that  explosion  Pn/Lg  indeed  is  higher  than 
earthquake  Pn/Lg  for  all  of  the  regionalized  Chinese  velocity  models,  we  should  keep  in  mind  that 
this  difference  is  much  smaller  than  the  variations  in  Pn/Lg  ratios  caused  by  propagation 
differences,  which  may  differ  by  as  much  as  a  factor  10.  Attempts  to  discriminate  between 
underground  nuclear  explosions  and  earthquakes  based  on  Pn/Lg  ratios  in  the  passband  of  our 
synthetics  should  therefore  be  undertaken  with  care.  For  the  broadband  frequency  range  the 
explosion/earthquake  quotient  can  be  both  smaller  or  larger  than  1  depending  on  the  velocity 
structure,  making  this  frequency  range  poorly  suited  for  discrimination  purposes. 

5.  Conclusions 

This  study  has  shown  that  the  PP  phase  can  effectively  be  used  to  model  crustal  and  upper 
mantle  P  wave  velocity  characteristics.  Long-period  PP  waveforms  can  be  clearly  discerned 
beyond  28°;  in  broadband  data  this  phase  can  be  identified  at  distances  down  to  25°.  The  strength 
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of  precursors  to  PP  can  be  used  to  interrogate  the  lid  velocity  gradients.  Long-period  direct  P 
waveforms  aie  not  as  sensitive  to  changes  in  upper  mantle  velocity  as  the  PP  phase,  but  details  of 
observed  P  triplication  behavior  in  broadband  data  can  be  used  to  put  constraints  on  the  velocity 
structure. 

Long-period  modeling  has  revealed  that  the  average  velocity  structure  of  the  crust  and  upper 
mantle  under  China  is  distinctive  in  character,  and  is  intermediate  between  models  for  active 
tectonic  regions  and  for  stable  shields.  Our  preferred  average  model,  WCH,  has  a  50  km  thick 
crust  with  a  velocity  of  6.4  km/s,  a  Pn  velocity  of  8.15  km/s  and  a  low  velocity  zone  between  100 
km  and  175  km  depth.  Velocities  at  greater  depths  are  similar  to  those  of  other  upper  mantle 
velocity  models.  This  model  can  match  the  majority  of  observed  long-period  waveforms  in  the 
14O-40O  distance  range  across  China.  However,  with  its  thick  crust  and  relatively  high  Pn  velocity 
it  is  more  characteristic  of  the  fold  belt  region  of  western  China,  which  is  preferentially  sampled  by 
most  of  the  paths  used  in  this  study.  These  two  features  contrast  with  the  thinner  crust  and  lower 
Pn  velocities  that  are  commonly  observed  for  active  tectonic  regions. 

Modeling  of  broadband  waveforms  gives  evidence  for  significant  lateral  deviations  from 
model  WCH.  Paths  traversing  eastern  China  indicate  a  constant  velocity  lid  of  about  100  km 
thickness,  an  observation  in  agreement  with  results  of  Holt  &  Wallace  (1990)  for  southeastern 
China.  The  Pn  velocity  in  this  region  ranges  from  8.12-8.14  km/s  and  crustal  thickness  from  27- 
35  km.  For  northern  China  crustal  thickness  is  between  43-50  km  and  a  Pn  velocity  of  about  8.0 
km/s.  There  is  some  preliminary  evidence  for  a  strong  lid  with  a  positive  velocity  gradient  of  about 
100  km  thickness.  For  the  Tibetan  plateau  we  observe  greater  crustal  thickness  (50-60  km)  and 
higher  Pn  velocity  (8.15-8.25  km/s)  in  combination  with  a  more  pronounced  low  velocity  zone 
than  that  of  the  average  model  WCH.  Lower  velocities  may  extend  below  200  km  depth.  These 
results  are  similar  to  those  of  the  S-wave  velocity  study  of  Lyon-Caen  (1986)  and  indicate  a  thin 
weak  lid  of  about  75  km  thickness.  We  infer  from  our  and  previous  results  that  crustal  shortening 
and  thickening  is  the  most  likely  process  to  have  given  the  Tibetan  plateau  its  thick  crust  and  high 
average  elevation.  The  Indian  plate  has  probably  only  slid  under  the  southernmost  part  of  the 
plateau.  Mantle  convection  associated  with  the  crustal  shortening  process  has  led  to  strong  lateral 
heterogeneity  in  the  upper  mantle  under  Tibet.  Higher  Pn  velocities  in  western  Tibet  confirm 
anomalously  high  Rayleigh  wave  phase  velocities  found  by  Brandon  &  Romanowicz  (1986)  and 
might  indicate  that  this  region  is  presently  undergoing  crustal  shortening. 

Variations  in  crustal  and  upper  mantle  velocity  structure  under  China  have  a  strong  influence 
on  the  propagation  of  the  Pn  and  Lg  phases;  Pn'/Lg  ratios  vary  by  as  much  as  an  order  of 
magnitude  for  the  regionalized  Chinese  velocity  models.  Of  the  model  variations  explored,  crustal 
thickness  is  the  most  important  factor  governing  the  propagation  of  Lg;  Pn  rms  amplitudes  are 
strongly  affected  by  crustal  thickness  and  variations  in  the  lid  velocity  gradient  Baseline  shifts  in 
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velocities  of  the  uppermost  mantle  have  only  a  minor  effect  on  the  propagation  of  Pn  and  Lg. 
Although  we  did  not  explore  the  effect  of  variations  in  crustal  velocity  this  is  also  likely  to  be  an 
important  factor  in  the  propagation  of  these  phases.  Short-period  Lg  has  very  stable  propagation 
characteristics:  high  frequency  Pn  rms  amplitudes  exhibit  a  more  complicated  behavior.  For  the 
regionalized  Chinese  velocity  models,  short-period  Pn/Lg  ratios  for  the  explosive  source  are 
consistently  higher  than  those  for  the  earthquake  source  However  differences  in  this  ratio  between 
the  explosion  source  and  the  earthquake  source  are  considerably  smaller  than  variations  in  Pn/Lg 
ratio  caused  by  propagation  differences.  Discrimination  attempts  based  on  this  ratio  should 
therefore  be  undertaken  with  care.  Broadband  Pn/Lg  ratios  for  the  explosion  source  can  be  both 
higher  or  lower  than  those  for  the  earthquake  source  making  this  passband  poorly  suited  for 
discrimination  purposes. 
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Table  1.  Earthquake  source  parameters 


event 

no 

date 

lat 

(deg) 

Ion 

(deg) 

depth 

(km) 

strike/dip/rake 

(degs) 

dur 

(sec) 

mb 

1 

07/23/88 

48.72 

90.51 

13 

345/72/1802 

2.8 

5.5 

2 

08/20/88 

26.78 

86.61 

57 

230/23/21 

3.0 

6.4 

3 

04/15/89 

29.98 

99.24 

10 

273/29/-731 

5.0 

6.2 

4 

04/20/89 

57.14 

121.92 

28 

6/56/1701 

4.7 

6.0 

5 

04/25/89 

30.04 

99.45 

10 

273/29/-732 

2.4 

6.1 

6 

05/03/89 

30.07 

99.48 

17 

273/29/-732 

5.2 

5.9 

7 

01/14/90 

37.77 

91.89 

8 

302/82/1052 

2.0 

6.1 

8 

03/05/90 

36.85 

73.01 

17 

192/36/-461 

5.0 

5.8 

9 

06/14/90 

47.89 

85.12 

32 

33/75/12 

6.5 

6.2 

10 

08/03/90 

47.95 

84.96 

■  34 

113/53/1731 

2.4 

6.1 

11 

01/05/91 

23.48 

95.98 

13 

20/70/1102 

3.0 

6.3 

12 

10/19/92 

30.85 

78.79 

16 

112/78/841 

7.0 

6.4 

13 

08/19/92 

42.10 

73.54 

30 

245/34/701 

6.0 

6.8 

Event  locations  and  magnitudes  from  NEIC,  depths  and  source  durations 
(dur)  from  waveform  modeling,  focal  mechanisms  from  1)  Harvard  CMT 
and  2)  waveform  modeling. 


Table  2.  Model  WCH 


depth(km)  vel(km/s) 

depth(km)  vel(km/s) 

0.000 

6.400 

270.000 

8.500 

50.000 

6.400 

280.000 

8.540 

50.000 

8.150 

290.000 

8.580 

100.000 

8.150 

300.000 

8.620 

110.000 

8.140 

310.000 

8.645 

120.000 

8.125 

320.000 

8.670 

130.000 

8.110 

330.000 

8.695 

140.000 

8.090 

340.000 

8.720 

150.000 

8.065 

350.000 

8.747 

160.000 

8.050 

360.000 

8.775 

170.000 

8.050 

370.000 

8.799 

180.000 

8.100 

380.000 

8.823 

190.000 

8.165 

390.000 

8.847 

200.000 

8.215 

400.000 

8.872 

210.000 

8.262 

405.000 

8.883 

220.000 

8.305 

420.000 

8.910 

230.000 

8.347 

420.000 

9.240 

240.000 

8.380 

431.000 

9.292 

250.000 

8.420 

442.000 

9.333 

260.000 

8.460 

452.000 

9.375 
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Table  3.  Value  of  b  for  fit  of  y=a-xb  to  Lg  decay  curves 


WCH 

east  China 

north  China 

-2.38 

-2.33 

-2.46 

-2.35 

broadband 

-1.70 

-1.68 

-1.08 

-1.36 

-2.41 

-2.42 

-2.09 

-1.81 

broadband 

-1.69 

-1.73 

-1.81 

-1.26 

Table  4.  Value  of  b  for  fit  of  yga-xt>  to  Pp  decay  curve 


source  type 

WCH 

east  China 

broadband 

-2.78 

-2.86 

-2.44 

broadband 

-2.80 

-2.79 

-2.41 
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Figure  1.  Tectonic  structure  sketch  of  China.  Outlines  of  tectonic  provinces  are  from  Mooney 
(1992).  Basin  and  platform  material  is  shown  in  grey  scale,  active  fold  belts  in  white.  The 
Chang  Thang  region  in  nonh-central  Tibet  which  is  characterized  by  Quaternary  volcanism  is 
indicated  by  the  hatched  lines. 
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Figure  2.  Map  showing  event  locations  (stars),  CDSN  broadband  stations  (circles  with  station 
codes)  and  paths  used  in  this  study.  Tectonic  areas  are  the  same  as  in  Figure  1.  Event  numbers 
correjpond  with  those  of  Table  1. 


velocity  (km/s) 


Figure  3.  Suite  of  P  wave  velocity  models  used  for  initial  long-period  modeling. 
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Figure  4.  Comparison  between  observed  and  synthetic  long-period  P  and  PP  waveforms  for  the 
velocity  models  shown  in  Figure  3  at  PP  triplication  distances.  Model  WCH  best  matches 
general  characteristics  of  these  waveforms  for  paths  traversing  different  tectonic  regions.  The 
first  arrival  of  PP  energy  as  picked  from  the  data  is  indicated  by  the  dashed  vertical  line.  The 
corresponding  paths  can  be  identified  from  Rgure  2. 
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PP  distance  (deg) 


P  distance  (deg) 

Figure  5.  P  and  PP  travel-time  curve  calculated  for  the  first  order  model  WCH  and  a  surface 
source.  The  reduction  velocity  is  10  km/s.  Triplication  branches  associated  with  the  420  and  670 
km  discontinuities  are  denoted  in  their  normal  convention.  The  dashed  branch  results  from  the 
low  velocity  zone  triplication  and  has  negligible  energy  associated  with  it 


Figure  6.  Comparison  between  observed  and  synthetic  long-period  P  triplication  waveforms  for 
the  velocity  models  shown  in  Figure  3  illustrating  that  model  WCH  can  also  match 
characteristics  of  these  waveforms  for  paths  throughout  die  China  region. 
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Figure  7.  Similar  comparisons  to  those  of  Figures  4  and  6  showing  that  deviations  from  the 
average  velocity  model  WCH  may  be  expected  in  the  upper  mantle  under  China.  The  first  arrival 
of  PP  energy  as  picked  from  the  data  is  indicated  by  the  dashed  vertical  line.  Branches  of  the  P 
triplication  arrives  correspond  to  those  in  Rgure  5. 


Figure  8.  Comparison  of  observed  and  synthetic  broadband  P  and  PP  waveforms  for  paths  mainly 
traversing  the  Tibetan  plateau.  This  illustrates  that  precursors  to  PP  can  be  used  to  infer  lid 
velocity  gradients.  Models  used  to  generate  synthetics  are  shown  in  the  velocity  profile. 


Figure  9.  Similar  comparison  as  Figure  8  for  broadband  P  waveforms  to  show  that  the  upper 
mantle  under  Tibet  can  be  characterized  by  higher  Pn  velocities  and  lower  velocities  below  200 
km  depth  relative  to  the  average  model  WCH. 


Figure  10.  Same  as  Figure  8  but  now  for  paths  traversing  eastern  China.  This  region  is 
characterized  by  a  constant  velocity  lid  and  Pn  velocities  that  are  somewhat  lower  than  Aose  of 
model  WCH. 
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Figure  11.  Similar  comparison  to  that  of  Figure  8  for  paths  traversing  the  northern  Chinese  and 
Mongolian  fold  belts.  These  paths  indicate  a  positive  velocity  gradient  in  the  lid  in  combination 
with  lower  Pn  velocities  than  those  of  model  WCH. 
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Figure  12.  Comparison  of  the  regionalized  P  wave  velocity  models  derived  in  this  study  with  P 
and  S  wave  velocity  models  from  previous  studies. 

34 


Figure  13.  Summary  of  crustal  and  upper  mantle  characteristics  under  China  at  4  depth  ranges  as 
derived  from  broadband  waveform  modeling.  Tectonic  areas  correspond  to  those  in  Figure  1. 
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Figure  13  cont. 
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reduced  time  (sec) 


Figure  14.  Synthetic  broadband  waveforms  generated  with  model  WCH  and  the  explosive  source. 
All  traces  are  normalized  to  their  maximum  amplitude.  The  reduction  velocity  of  8.2  km/s  aligns 
the  Pn  arrivals  at  almost  constant  time.  The  Pn  and  Lg  windows  used  for  the  rms  amplitude 
calculations  are  bounded  by  the  dashed  lines. 
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Figure  15.  Spectra  of  the  broadband  synthetic  waveform  at  7°  of  Figure  14,  and  the  same 
waveform  filtered  with  the  OBN  Kimos  insttumenL  This  illustrates  the  effea  of  the  short-period 
filtering  on  the  frequency  content  of  the  synthetic  waveforms. 
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Figure  16.  Same  as  Figure  14  but  now  for  the  waveforms  filtered  with  the  short-period  response. 
Long-period  Rayleigh  wave  modes  have  dissappeared  from  the  Lg  window  and  the  Pn  arrival 
has  become  more  impulsive  and  complex. 
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distance  (km) 

Figure  17.  Short-period  rms  Pn  and  Lg  amplitudes  and  their  ratios  versus  distance  for  4  different 
P  wave  velocity  models  and  the  explosive  source.  The  velocity  models  can  be  identified  in 
Figure  12;  the  synthetic  waveforms  for  model  WCH  are  shown  in  Figure  16. 
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Figure  19.  Same  as  Figure  17  but  now  for  die  earthquake  source. 
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Figure  21.  Quotient  of  explosion  and  earthquake  Pn/Lg  ratios.  Though  the  short-period  quotient  is 
consistently  higher  than  1  for  the  regionalized  Chinese  velocity  models  it  is  considerably  smaller 
than  Pn/Lg  variations  due  to  propagation  differences  as  can  be  inferred  from  Figures  17-20. 
Depending  on  the  velocity  structure,  the  broadband  quotient  can  be  both  higher  and  lower  than 
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Appendix  A 


Overview  of  long-period  modeling  results.  Data  are  shown  in  top  trace  and  are  underlain  by 
synthetics  for  models  S25,  K8,  WCH  and  T7.  The  synthetic  for  the  preferred  model  is  shown 
solid,  other  synthetics  are  dashed.  Data  and  synthetic  are  aligned  on  their  first  motions  and  are 
normalized  relative  to  their  peak  amplitude.  Event  numbers  and  stations  correspond  to  those  in 
Figure  2  and  Table  1. 
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Appendix  B 

Overview  of  broadband  modeling  results.  Top  traces  show  match  between  the  observed  waveform 
and  the  synthetic  generated  for  model  WCH,  the  bottom  trace  gives  the  best  match  obtained  for  this 
waveform.  The  velocity  profiles  indicate  changes  made  to  model  WCH  to  obtain  this  match.  Event 
numbers  and  stations  correspond  to  those  in  Figure  2  and  Table  1. 
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Appendix  C 

Overview  of  synthetic  regional  waveforms  generated  to  explore  the  effect  of  the  regionalized 
Chinese  velocity  models  (Figure  12)  on  the  propagation  of  the  Pn  and  Lg  phases.  Synthetics  are 

calculated  in  the  10-10°  distance  range  at  1°  intervals.  The  waveforms  are  plotted  with  a  reduction 
velocity  of  8.2  km/s  and  are  normali^  to  their  peak  amplitudes. 
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